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Abstract 

The dielectric behavior of a linear cluster of two or more living cells connected by tight junctions 
is analyzed using a spectral method. The polar izability of this system is obtained as an expansion 
over the eigenmodes of the linear response operator, showing a clear separation of geometry from 
electric parameters. The eigenmode with the second largest eigenvalue dominates the expansion as 
the junction between particles tightens, but only when the applied field is aligned with the cluster 
axis. This effect explains a distinct low-frequency relaxation observed in the impedance spectrum 
of a suspension of linear clusters. 
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I. INTRODUCTION 



Particle polarizability governs the electric response for many inhomogeneous systems 
ranging from biological cells to plasmonic nanoparticles and depends strongly on both its 
dielectric and geometric properties. Analytical models have been reported l[ Q only for 
spherical and ellipsoidal geometries, whereas more complex geometries have been approached 
by direct numerical solution of the field equations using, for example, the finite difference 
methods [3], the finite element method 4|, the boundary element method 

aa 

, or the 

boundary integral equation (BIE) {7]. 

In a simplified representation, biological cells can be regarded as homogeneous particles 
(cores) covered by thin membranes (shells) of contrasting electric conductivities and permit- 
tivities. Complex geometries occur when cells are undergoing division cycles (e.g. budding 
yeasts) or are coupled in functional tissues (e.g. lining epithelia or myocardial syncytia). 
In these cases, the dielectric/impedance analysis of cellular systems is far more complicated 
than previous models 8-1^, which considered suspensions of spherical particles. Intriguing 



dielectric spectra 



llj reveal distinct dielectric dispersions with time evolutions consistently 



related to tissue functioning or alteration, identifying a possible role of cell connectors (gap 
junctions) in shaping the overall dielectric response. 

A direct relation between the microscopic parameters and experimental data can be 
analytically derived only for dilute suspensions of particles of simple shapes, and is rather 
challenging for system with more realistic shapes, where only purely numerical solutions have 
been available. In this work we demonstrate that a spectral representation of a BIE provides 
the analytical structure for the polarizability of particles with a wide range of shapes and 
structures. The numerically calculated parameters encode particle's geometry information 
and are accessible by experiments. 

By using single and double-layer potentials |l2| , the Laplace equation for the fields inside 
and outside the particle is transformed into an integral equation. A spectral representation 
for the solution of this equation is obtained providing the eigenvalue problem for the linear 
response operator is solved. Although not symmetric, this operator has a real spectrum 
bounded by -1/2 and 1/2 13Nl5| and its eigenvectors are orthogonal to those of the conjugate 
double-layer operator. A matrix representation is obtained by using a finite basis of surface 
functions. 



2 



The true advantage of the spectral method is that the eigenvalues and eigenvectors of the 
integral operator provide valuable insight into the dielectric behavior of clusters of biological 
cells. The eigenvectors are a measure of surface charge distributions due to a field. Only 
eigenvectors with a non-zero dipole moment contribute to the polarizability of the particle. 
We call these dipole-active eigenmodes. An effective separation of the geometric and mor- 
phologic properties from dielectric properties is therefore achieved 16| . We also show that 
for a particle covered by multiple confocal shells, the relaxation spectrum is a sum of Debye 
terms with the number of relaxations equal to the number of interfaces times the number 
of dipole-active eigenvalues. This is a generalization of a previous result [17| on cells of 
arbitrary shape. 

Our method is related to another spectral approach which uses an eigenvalue differentia 
equation 



1 3M2 if] . This method has been applied to biological problems by Lei et al. 



22] 



and by Huang et al. 23J. These authors, however, considered homogeneous cells with much 
simpler expression for cell polarizability. The BIE spectral method seeks a solution on the 
boundary surface defining the particle, as opposed to the eigenvalue differential equation, 
where the solution is defined in the entire space. 

In a previous study on double (budding) cells it was shown that before cells separation an 
additional dispersion occurs 24]. Moreover, in recent papers {3, 25] numerical experiments 
have shown that the dielectric spectra of a suspension of dimer cells connected by tight 
junctions exhibit an additional, distinct low-frequency relaxation. Our numerical calculation 
shows that the largest dipole-active eigenvalue approaches the value of 1/2 as the junction 
become tighter. Although the coupling of this eigenmode with the electric field stimulus 
is relatively modest (the coupling weight is about 1-2 %), this eigenmode has a significant 
contribution to the polarizability of clusters. Thus the eigenmodes close to 1/2 induce an 
additional low-frequency relaxation in the dielectric spectra of clustered biological cells even 
though the coupling is quite small. Needle-like objects, such as elongated spheroids or long 
cylinders, have similar polarizability features. 

In this paper we consider rotationally symmetric linear clusters made of up to 4 identi- 
cal particles covered by thin insulating membranes and connected by junctions of variable 
tightness. Convenient and flexible representations for the surfaces describing these objects 
are provided. The number of relaxations in the dielectric spectrum of the linear clusters, 
their time constants and their relative strengths are analyzed in terms of the eigenmodes of 
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the linear response operator specific to the given shape. 



II. THEORY 

A. Effective permittivity of a suspension 

We consider a suspension of identical, randomly oriented particles of arbitrary shape and 
dielectric permittivity, ei, immersed in a dielectric medium of dielectric permittivity, Eq. 
The dielectric permittivities are in general complex quantities and the theory described here 
applies also for time-dependent fields, providing that the size of a particle is much smaller 
than the wavelength. When an applied uniform electric field interacts with the suspension, 



the response of the system is linear with the applied fie 



the whole sample can be measured and is defined by [25 



d and an effective permittivity for 



27J: 



£sus — T J _ j, ^ ■ \ L ) 
J 3 

This result is obtained in the limit of low concentration, weak intensity of the stimulus field, 
and using an effective medium theory within the dipole approximation. Here / = NVi/V is 
the volume fraction of all N particles, each of volume Vi, with respect to the total volume 
of the suspension V. The averaged normalized polarizability a of a particle is defined as 
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a = T^7 / (^^)B(N)-NdQ N dV (2) 
47rVi J Vl Jn N \ e J 

where E (N) is the electric field perturbation created inside the particle under a normalized 
applied electric excitation with direction N and cIQn is the solid angle element generated 
by that direction. The above normalized polarizability is dimensionless and is obtained by 
multiplying the standard polarizability of a particle with the factor Air/V\. In the following 
we will refer only to normalized polarizability, thus, without any confusion, the normalized 
polarizability a will be simply called polarizability. The directional average in ([2]) is equiv- 
alent to the averaged sum over three orthogonal axes due to the fact that the problem is 
linear with respect to the applied field. The latter is more convenient from computational 
point of view. 

The electric field inside a particle is obtained by solving the following Laplace equation 
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for the electric potential $: 



A$ (x) =0, x G K 3 \S 

$| + = $]_ , x G S 

£ off | + = £ iff |_ > x e S 

$ — > —x • N, Ixl — > oo 



(3) 



where !ft 3 is the euclidian 3-dimensional space and £ is the surface of the particle. The 
derivatives are taken with respect to the normal vector n to the surface S. 

Due to the mismatch between the polarization inside and outside the object, electric 
charges accumulate at the interface £ and create an electric potential which counteracts the 
uniform electric field stimulus. The solution of the above Laplace problem ([3]) is therefore 
formally given by 

d>(x) = -x - N + J- / dS(y) (4) 

47r Jt, l x - y| 

The single layer charge distribution fi induced by the normalized electric field is a solution 
of the following BIE, obtained by inserting solution (TjJ in equations (J3]) 

/i(x) If n(x) • (x - y) 



2A 4vr 



/ MY) ' (X ,3 y) dE(y) = ^(x) • N (5) 
Jt, l x — y| 



Here the parameter A = (^i — £q) / (si+eo) isolates all the information regarding the dielectric 
properties for this problem. 

On using the linear response operator M that acts on the Hilbert space of integrable 
functions on the surface E, 

MM = i / My) n( f ' (x " y) dE(y), (6) 

the integral equation (j^J) is written as 

(l/(2A)-M)/i = n-N. (7) 

The integral operator ([6]) is the electric field generated by the single layer charge distribution 
\i along the normal to the surface. It encodes the geometric information and has several 



interesting properties 
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151 ]. Its spectrum is discrete and it is not difficult to show that 
all of its eigenvalues are bounded by the [-1/2, 1/2] interval. Although non-symmetric, the 
operator ([6]) has real non-degenerate eigenvalues. The eigenvectors are biorthogonal, i. e., 
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they are not orthogonal among themselves, but orthogonal to the eigenvectors of the adjoint 
operator 

MVH^jU) "'^VW ), (8) 

which is associated with the electric field generated by a surface distribution of electric 
dipoles (double layer charge distribution). Therefore, if \uu) is a right eigenvector of M 
corresponding to eigenvalue Xk, M\uk) = Xk\uk) and (tv| is a left eigenvector corresponding 
to eigenvalue Xk', {vk'W^ = Xk'(vk'\, then 

(vk'\uk) = 5k>k, (9) 

with the scalar product defined as the integral over the interface S, 

(/i|/ 2 >= / /r(x)/ 2 (x)dS(x). (10) 

The value 1/2 is always the largest eigenvalue of the operator M, regardless the geometry 
of the object. This is immediately seen if the object is considered to be conductor {e\ — >■ oo), 
and then the interior electric field has to be zero. In that case A = 1, and the charge density 
that generates a vanishing internal electric field obeys the equation (1/2 — M)/i = 0, and 
therefore 1/2 is an eigenvalue of M. However, this eigenmode is not dipole-active and does 
not contribute to the total polarization of the object. The operator ([6]) is insensitive to a 
scale transformation, which means that its eigenvalue and eigenvectors depend only on the 
shape of the object and not on its size, or electrical properties. 

By employing the spectral representation of the resolvent of the operator M 

{z - M)- 1 = J> - XkY X \u k ){v h \, (11) 

k 

the solution of equation (J7|) is obtained for z = 1/(2A) as 

ii = 5> fe |n ■ N)(1/(2A) - XkVW). (12) 



The polarizability of the homogeneous particle is obtained by using the distribution ( fl2l) 
to build the solution (Hj) of the Laplace equation and use it in equation (j2J). It has been 
shown that, operationally, the polarizability is simply the dipole moment of the distribution 



(fl2"j) over unit volume [25|, [27] 
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1 1 ^ (x • Ni\u k ) {v k \n - Ni) 
"=3^5 1/(2A) _ X , • 

where Nj are three mutually orthogonal vectors (directions) of unit norm. The factor 
(1/(2A) — Xk)^ 1 is a generalized Clausius-Mosotti factor. Each dipole-active eigenmode 
contributes to a according to its weight p k = |-^Xli( x ' (vk\n ■ Nj), which deter- 

mines the strength of coupling between the uniform electric field and the fc-th eigenmode 
and contains three components Pk,i = (x • Ni|iXfc) (vk\n ■ Ni)/Vi. Equation (TT5]) shows a 
clear separation of the electric properties, which are included only in A, from the geometric 
properties expressed by Xk and p k - 

B. Shelled particles 

The polarizability of an object covered by a thin shell with permittivity e$ can be calcu- 
lated in a similar fashion. The electric field is now generated by two single layer distributions, 
and boundary conditions are imposed twice, for Si and for E 2 . The surface is the outer 
surface of the shell and E 2 is the interface between the particle and the shell. The solution 



of a shelled particle in terms of single layer potentials has the form 28| 



*(x) = -x ■ N + -L / <E(y) + -L f {M. ffibr), (14) 

4tt y Sl |x-y| 4vr J^ 2 |x-y| 

where \i\ and /z 2 are the densities defined on surface Ei and E 2 , respectively. Four integral 
operators M u , M 12 , M 2 \ and M 22 are defined, depending on which surface are variables x 
and y. For example Mu is defined when x and y are both on Si, M\% is defined by x on 
Si and y on E 2 , and so on. Thus 

MM = ^zf h(y) n(x l ' ( * ,7 y) rfs (y) ( 15 ) 



4vr i s . |x - y 

for i, j = 1, 2. The equations obeyed by /xi and /i 2 are 

/ii/(2Ai) -Mnf/ii] -M 12 [// 2 ] = n-N 

(16) 

^ 2 /(2A 2 ) - M 2 i[^i] - M 22 [/i 2 ] = n • N. 
Here the electric parameters are: Ai = (eg — £ )/(£g + £0)5 and A 2 = (ei — £g)/(ei + £5). 

We further assume a confocal geometry, i. e. the surface Ei is a slightly scaled version 
of E 2 , with a scaling factor 77 close to unity. This assumption does not provide constant 



thickness for the shell, but our main results should remain at least qualitatively valid 



In the 

can show 



imit of very thin shells, and using the scaling properties of the operator M, one 



3,0 



that all four M operators are related to M — Mn 
MvM =-n-\n/2 + M\n\) 

M u M = -/i/2 + MM (17) 
M 22 [fi) =M\pi\. 

Equations (FIB]) can then be arranged in a matrix form as 



l/(2Ai-M) (l/2 + M)/r/ 3 
-1/2 + M 1/(2A 2 )-M 



)l 






( n-N\ 








l n-N J 



(18) 



By knowing the eigenvectors and the eigenvalues of M the charge densities fi\ and [i 2 can 
be found by inverting the matrix in (fl8|) . For example, /ii is 

^(ifc - Xk) + (| + Xk) 

/^Z^ - ^ \ n 77^ -(«fe|n-N)|« fc ). (19) 

* ^ [m ~ Xk ) \2x-2 ~ + (2 + Xk) (2 - Xfc) 

The field generated by the the two distributions H\ and fi 2 outside the particle is the 
same as the field generated by an equivalent single layer distribution 



He 



£>k|n ■ N)(l/(2A fc ) - X/O^K), (20) 



k 



where Xk = (£&— £0) / (£k+£o) and the equivalent permittivity 6k is defined for each eigenmode 



as: 

£1 ~ £s 



ik £S \ 1 + e s + 5(1/2 - Xk )e l + 6(1/2 + X k)e s )> (21) 
where 5 = rf — 1 <C 1. The distribution ( 120]) is similar with the distribution f JT2|) obtained for 
a homogeneous particle, except that A has to be replaced for each mode with an equivalent 
quantity \\. Equation (}2"T]) can be applied recursively for a multi-shelled structure. The 
strict separation of electric and geometric properties is weakened in this case, because the 
shape-dependent eigenvalue Xk appears now in the electric equivalent quantity X k . 

The polarizability of the shelled particle is obtained by using the distribution ( 1201 to 
build the solution (j3J) of the Laplace equation and use it in equation (JSJ, to get 

_ 1 1 y (x' NjjtXfc) (v k \n ■ Nj) 
^'^Yk 1/(2A*)-X* 
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Equation ( 1221) is obtained by replacing A with Afc in Eq. ( JT3l) . The parameter Vi in Eq. ( 1221) 
is the total volume of the cell (the core and the shells). In the limit of a dilute suspension 
of identical shelled particles, with a low volume fraction /, the effective permittivity (TjQ) is 

V k ( 1 / 2 + Xk)eo + {l/2-Xk)e k J 

C. The Debye relaxation expansion 

In general, the effective permittivity e sus of a suspension of objects with m shells will have 
m+1 Debye relaxation terms for each dipole active eigenmode. The proof is recursive and is 
based on partial fraction expansion with respect to variable iu of equations ( 122]) and (123]) . 
provided that the complex permittivity of various dielectric phases is e = e — ia/(ue vac ) 
where i = a/— 1 and e vac is the permittivity of the free space (8.85 x 1CT 12 F/m). Thus the 
first Debye term comes out from (1231) and the remaining m Debye terms result from (T2TJ 
by the homogenization process described for shelled particles. Hence, a suspension of cells 
with m shells (and m + 1 interfaces) has a dielectric spectrum containing a number of Debye 
terms equal to m + 1 times the number of dipole-active eigenvalues. 

The suspension effective permittivity e sus has the expansion 

e sus = e f + A£ kj/(l + iwT kj ) (24) 

where ef = £hf — icrif/(we mc ), £m is the high-frequency permittivity, and an is the low- 
frequency conductivity; Aejy and are the dielectric decrement and the relaxation time 
of the kj Debye term, respectively; index k enumerates the dipole-active eigenmodes and 
index j enumerates interfaces. 

Although the measurable bulk quantities in equation (1241) are directly correlated with the 
microscopic (electric and shape) parameters, a solution of the inverse problem, which aims 
at obtaining the microscopic information non-intrusively, from the effective permittivity, is 
in general difficult, if not impossible for the general multi-shell structure. However, biolog- 
ical cell has a thin and almost non-conductive membrane, and several simplifications and 
approximations can be made. Two Debye relaxation terms in the effective permittivity e sus 
are expected for each dipole-active eigenmode, corresponding to the two interfaces which 
define the membrane. 
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The first relaxation is derived from the equivalent permittivity ([21]) which can be written 
also as Debye relaxation terms: 



h = e + Ae/(l + iuT). (25) 

The relaxation time T that is given by the poles of Ik in fl2T|) is a quite good approximation 
of the first relaxation time Tj.% 

rp _ (I + 5/2 + 5 X k)e s + 5 (1/2 - Xk )e 1 „ T (9f ,, 
~ Evae {l + 6/2 + 6 X k)<Ts + 6(1/2 -xk)<ri~ hl ' [ ' 

The main reason is as follows. At frequencies close to 1/T there is a huge change in Ik of 

order es/ (6(1/2 — Xk)), and consequently a significant change in the total permittivity e sus 

given by ( )23l) . Therefore T provides an approximate value for the relaxation time Tki of the 

suspension effective permittivity e sus . 

For a non-conductive shell as ~ 0, or more precisely when cr s <^6 (1/2 — Xk) the 

relaxation time ( ]26]) is: 

Tki « e««c es/(5 • <ri(l/2 - xfc)) (27) 

showing a strong dependence on the thickness of the shell and on the shape of the particle, 
through the eigenvalue Xk- Due to the small parameter 6 in ( 1271) the first relaxation (i. e., 
membrane relaxation) tends to have a lower frequency than the second relaxation, which is 
present even for particles with no shell (see the discussion below). In addition, cumbersome 
but straightforward calculations provide the dielectric decrement Ae^ in (|24|) 

Ae k i « Afp k e s (6(l/2 + Xfc ) 2 (l/2 - Xk))' 1 (28) 

that is very large due to the same strong dependence on the thickness of the shell. The 
effect is even more dramatic when the second largest eigenvalue is very close to the largest 
eigenvalue, (1/2 — % 2 ) — > 0, like in the case of two cells connected by tight junctions. 

For a suspension of shelled spheres rj = 1 + AR/R and 6 = t] 3 — 1 ~ 3AR/R, where 
AR is the thickness of the membrane, R is inner radius, and R + AR is the total radius. 
Thus both and Aeki are proportional to R and es and inverse proportional to AR like 
in the Pauly-Schwan theory [l|, Q]- Moreover, the dielectric decrement Aski in (128]) is a 
generalization of equation (54a) in Ref. 8J. In the same time, the relaxation time (]27j) 
differs with respect to equation (56a) in Ref. js] only by the conductivity term. We will 
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show elsewhere that a more appropriate treatment of the relaxation times recovers also the 



relaxation time given by equation (54a) in Ref. 

Thus, a non-conductive and thin shell/membrane produces a large relaxation of the com- 



plex permittivity of the suspension 3l|. The experimental evidences further support these 
theoretical facts: when attacking the membrane with a membrane disrupting compound ( for 
example a detergent) the relaxation almost vanishes as the cellular membrane is permeated 



32] 



For frequencies higher than l/X^ the cell permittivity is essentially determined by the 
dielectric properties of the cytoplasm, and does not depend on membrane's properties. The 
second Debye relaxation occurs at higher frequencies than the first (membrane) relaxation, 
and has the relaxation time 

rp „ (l/2 + x*)eo + (l/2-Xfc)gi 

I ^-^ {1/2 + Xk)ao + {1/2 _ Xk)ai W 
derived from the pole of equation ( 123|) . The corresponding dielectric decrement is 

Ae fc2 ~ /p*(l/2 - Xk)(£i(?o ~ £o^i) 2 x 
((l/2 + Xfe ) £o + (l/2-x^i) _1 x 
((l/2 + Xfc )^o + (l/2-X fc )a 1 )- 2 

The last two equations are similar to the ones that are given for spherical particles in Ref. 
s| (equations (46) and (49) in the aforementioned reference). The relaxation given by 
is basically the relaxation of a homogenous particle embedded in a dielectric environment 
and was also discussed in Ref. [22| by a closely related spectral method. If the conductivity 
of the cytoplasm is comparable to the conductivity of the outer medium, the decrement of 
the second relaxation is small such that it cannot be distinguished in the spectrum. On the 
contrary, if the conductivity of the outer medium is much greater or smaller that that of 
cytoplasm, than a second observable relaxation occurs. Unlike the membrane relaxation, 
this second relaxation depends only weakly on the shape. By assuming that <jq <C <j\ 
and by using a finite-difference method, this resonance was also obtained in [3lj and it was 
instrumental in explaining the experimental data on the fission of yeast cells of Asami et al. 



[33| by Lei et al. [22 1 . 



The shape of the particle is important because it affects the number of dipole-active 
eigenvalues and their strengths. In principle, each dipole-active eigenvalue introduces a new 
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relaxation in the dielectric spectrum, providing this relaxation is well separated from the 
others. A cluster with complex geometry can have several dipole-active eigenvalues, but 
unless the cluster is larger in one dimension then in the others, or there are tight junctions, 
the relaxations overlap to create broad features in the spectrum. An extra relaxation is 
introduced when the particles are covered by thin membranes. In addition, if (1/2 — Xk) 
for that eigenvalue, then the shell induced relaxation has low frequency, large relaxation 
time and large dielectric decrement. Based on the spectral BIE method, it is therefore 
possible to explicitly relate the dielectric spectra of cell suspensions to cell's geometry and 
electric parameters, and, even design fitting procedures to evaluate these parameters from 
measurements. 



III. RESULTS 

A. Numerical procedure 

The calculation of the effective permittivity for a suspension uses equations ([I]) or (1231) . 
and reduces then to finding the eigenvalues \k and eigenvectors \uk) and \vj.) of the linear 
response operator M. This problem is solved by employing a finite basis of NB functions 
defined on the surface S. A natural basis for a surface not far from a sphere is the generalized 
hyperspherical harmonics functions 

F im (x) = -±= Y lm (0(x), p(x)) , (30) 
V s ( x ) 

where s(x) is related to the surface element through dT, = s(x) dfl^ and dQ^ is the solid 
angle element. 



Another choice could be based on Chebyshev polynomials of the first kind 29] 



r lm (x) = -±= Ti (0(x)) e^W. (31) 
V s ( x ) 

Both bases are complete and orthogonal in the Hilbert space of square integrable functions 
defined on E. 

In this paper, we model the linear cluster of particles as an object with axial symmetry. 
We seek to find a surface of revolution for which the thickness of the interparticle joints can 
be varied without perturbing the overall shape of the object. We use two representations 
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for the surface £: (A) for clusters of two particles we use spherical coordinates {x,y,z} = 
{r(6) sin 9 cos </>, r(9) sin 9 sin <ft, r(9) cos 9}, and (B) for clusters with more than two particles 
we specify the surface in terms of a function g(z) as {x, y, z} = {g(z) cos0, g(z) sin0, z}. 
In the case B, the surface element is 



dH = g(z)\fl + g' 2 {z) dzdtp, 



(32) 



and the normal to surface £ is 



n 



l + g' 2 (z) 



^ cosy) ' 
sin if 

\ -9'{z) ) 



(33) 



In the basis of generalized hyperspherical harmonics the operator M has matrix elements 
given by 



2tt 



Mi 



lm;l'm' 



where 



and 



A(z, z' 



A(z, z',<p-<p')P[ n (cos9(z))P[r'(cos9(z')) j m &~^ G(z,z') dz dz' d<p dip' 

(34) 

G(z, z') = sjg{z)g{z')^/{r+g'{z)){l+g'{z')f 1 \ (35) 
{g{z) - g{z')) cos <p-(z- z')g'(z) 



(36) 



\g\z) + g\z<) - 2g(z)g(z>) cos0 + (z - z'ff 2 
After the angle integration in equation ( |34|) and by using the elliptic integrals given in 
the Appendix, the matrix elements are obtained by numerical evaluation of the resulting (z, 
z') double integral using an iVQ-point Gauss-Legendre quadrature [29|, |30(. Because of the 



integrable singularity apparent in the kernel of the operator M in equation ([6]), the mesh 
of z must be shifted from the mesh of z' by a transformation which insures that there is no 
overlap between the two meshes. 

The delta symbols S mm ' in equation (134"1) reflects the fact that we consider only objects 
with rotational symmetry in this paper. Moreover, for fields parallel with the cluster axis 
m = 0, while m = 1 for perpendicular fields. 

The convergence of the results is obtained in two steps. First, the number NQ of quadra- 
ture points is increased until the matrix elements of M converge, and then the size NB of 
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the basis set is increased until the relevant eigenvalues Xk and their corresponding weights 
Pk have acquired the desired accuracy. A necessary test for convergence is the fulfillment of 
the sum rules ^2 k Pk,i = 1 and ^ k XkPk,i = 1/2 with sufficient accuracy [18|, [l9|. Usually 
the convergence is fast in both the number of quadrature points and the size of basis, unless 
the system has a tight junction where some care must be considered in order to achieve the 
required accuracy of the eigenmodes with the eigenvalues close to 1/2. For a sphere there is 
just one dipole-active eigenmode which has eigenvalue \ = 1/6 and weight p — 1, while for 
an ellipsoid there is one dipole-active eigenmode along each axis. Fast and accurate solutions 
are achieved for spheroids with a basis size of NB =20 and with NQ = 64 quadrature points. 
In general, the size of the basis and the number of the quadrature points increase with the 
number of cells in the cluster and with the decreasing of the junction size. Thus, for our 
numerical examples a basis with iV5=35-40 and NQ =128 quadrature points are enough 
for a converged solution in the case of the dimers and NB=50 and NQ = 200 quadrature 
points in the case of the clusters with up to four cells. 



B. Two cells joined by tight junction 

The equation r(6) = (/i+cos 2 0)/(l — a cos 2 6) describes the shape of a two-particle cluster. 
Parameter h controls the tightness of the inter-particle junction and parameter a measures 
the deviation from a spherical shape. More precisely, h is the radius of the smallest circle 
at around the thinnest part of the junction. 

Figure [1] shows the effective permittivity for a suspension of particles with the following 
parameters: S\ = 70, <j\ = 0.25 S/m, es = 6, as = 0, Sq = 81, cr = 0.374 S/m, volume 
fraction / = 0.05, membrane thickness 5 = 0.00947275, and a = 0.2. The effective per- 
mittivity does not depend on the thickness parameter h when the stimulus electric field is 

uster axis. However, a new relaxation becomes apparent as h — > 0, 
|3|,|25|. 

Figure [2] presents the first 7 eigenvalues Xk and their weights P k ,i for a field parallel 
to the z axis. As the junction become tighter (h — > 0) more eigenvalues become dipole- 
active. While all eigenvalues are important in shaping the dielectric spectrum, the second 
largest eigenvalue X2 is crucial to explaining the occurrence of an additional relaxation at 



perpendicular to the c 
for parallel fields 



low frequencies, as observed for small h in |3j, |25( . Although its weight P2.1 also decreases for 



14 




FIG. 1: (Color online) The spectrum of the effective permittivity of a suspension of dimers with 
various junction thickness h, and with parallel and perpendicular field configurations. The sus- 
pension permittivity for an electric field perpendicular to cluster axis does not depend on h and is 
pointed by an arrow. 

small h, this dipole active eigenmode approaches 1/2 as the junction becomes tighter. Thus, 
according to equations ( |27|) and (1281) the effect of \2 is "enhanced" due to the presence of a 
nonconductive shell (as the case for biological particles analyzed in {3, 25 ] ) . Moreover, the 
decrease of P2,i is compensated by the increase of 1/(1/2 — x 2 ). 

The presence of a new relaxation at low frequency alon g w ith its relationship with the 
size of h has bee,, a.read y sing.ed out in Gheorghiu e t al. by us,ng the same method 
but without the analysis of dipole-active eigenmodes. Using a finite discrete model |3j, 
ihe relaxation was observed before the segregation during cell division, while other papers 
34, 35] fail to relate the size of h to the new relaxation, even though one of them jsjj 
employs essentially the same method as the one outlined in the present work. 

Figure |3] shows the charge distribution associated with the first four eigenvalues for two 
distinct values of h. The second eigenmode is an antisymmetric combination of net charge 
distributions (monopoles) on each particle of the dimer. The third charge distribution is an 
antisymmetric combination of charge distributions with a dipole moment on each part of the 
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FIG. 2: (Color online) The largest seven eigenvalues and their weights for a binary cluster with 
parameter a = 0.2, as a function of h. The inset shows the shape of the dimer. 

dimer and the forth distribution is antisymmetric combination of charge distributions with 
a quadrupole moment on each particle. At small h (tight junctions), charge accumulates in 
the vicinity of the junction 
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cos(e) 

FIG. 3: The first four eigenvectors for a dimer given by equation r{9) ; a = 0.2 and h = 0.5 (dotted 
line) and h = 0.01 (solid line). 

C. Clusters of more than two particles 

Smoother yet tight junctions would bring (1/2 — X2) closer to than sharp and tight junc- 
tions. The reason is simple: smoother junctions have the two parts of the dimer farther apart. 
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We have analyzed linear clusters of cells connected by smooth and tight junctions by using 
a (z, 4>) parameterization which describes a surface by {x = g(z) cos ip, y = g(z) simp, z). 
The construction starts from a dimer shape that resembles the shape of the epithelial cells 
like MDCK (Madin-Darby Canine Kidney) cells. An example of such shape, displayed in 
figure IU extends from — z max to z max and it can be decomposed in three parts: the left cap 
(-z mai < z < —Zi), the central part (— Z\ < z < Zi), and the right cap (z\ < z < z TaBX ). 
At position ±z± the shape function has its maximum. An m-cell linear cluster is ob- 
tained by repeating the central part m — 1 times and it extends from — L m to L m , where 
L m = z max + (m — 2)z±. Mathematically, the shape is described by: 

g(z + (m - 2)zt), for - L m < z < -L m + z max 
9m{z) = < g (mod(z + (1 + (-l) m )zi/2, 2z t ) - z t ) , for - L m + z max < z < L m - z max 
g(z - (m-2)z 1 ), for L m < z < L m - z max , 

(37) 

where mod(x, y) is the remainder of the division of x by y. For the examples considered 
here, the dimer shape function is: 

g(z) = 0.01 + 2.32317 z 2 - 11.9862 z 4 + 40.4045 z & 
-74.2226 z 8 + 79.142 z 10 - 51.8929 z 12 + 21.3096 z 1A (38) 
-5.35113 z 16 + 0.752147 z 18 - 0.045375 z 20 , 
with z max = 1.77377 and z x = z m ^/2. 

Tables [J and [Til list the most representative dipole-active eigenmodes for a trimer in 
perpendicular and parallel fields. Only the parallel field configuration has a dipole-active 
eigenvalue close to 1/2, with a relatively small weight. 

The results for linear clusters of up to four particles are displayed in Figure [5j The elec- 
tric parameters are the same as ones used for dimers in the previous section. An additional, 
distinct low-frequency relaxation emerges for clusters with more then one particles, only 
when the stimulus field is parallel with the symmetry axis. The relaxation frequency de- 
creases, while the intensity of these relaxations increases, as the number of cluster members 
increases. This behavior is explained again by the combination of eigenvalues close to 1/2, 
with thin non-conductive layers covering the cluster and is consistent with experimental data 
on ischemic tissues [ill ], which reports that the cell separation (closure of gap-junctions) is 
responsible for decrease and eventual disappearance of the low-frequency dispersion. 
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TABLE I: Most representative dipole-active eigenmodes and their weights for the trimer in parallel 
field. 



Xk 


Pk,i 


0.4996 


0.01305 


0.40642 


0.07769 


0.40448 


0.1391 


0.37763 


0.17366 


0.26409 


0.20519 


0.23809 


0.01839 


0.17557 


0.05104 


0.13522 


0.05234 


0.12165 


0.01066 


0.07192 


0.06197 


0.06649 


0.01239 


0.03479 


0.03899 


0.03362 


0.02823 


0.0281 


0.04688 


0.02377 


0.02519 



TABLE II: Most representative dipole-active eigenmodes and their weights for the trimer in per- 
pendicular field. 



Xk 


Pk,2 


0.1831 


0.64102 


0.06964 


0.01464 


0.05492 


0.01343 


0.0272 


0.01065 


0.01674 


0.05687 


0.01479 


0.13179 


0.00395 


0.05029 


-0.01635 


0.03246 
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FIG. 4: Smooth construction of a cluster (lower panel) from a dimer (upper panel). The parts 
determined by z £ [—Zx,Zi] are "glued" together with the ends of the dimer. The arrows show 
where the junctions will be placed in the cluster. 

In Figure [6] we plot Pfc ) i/(l/2 — Xk) versus (1/2 — Xk), which shows that the number of 
dipole-active eigenmodes increases with the number of particles in the cluster. According 
to (|27p and (|28|) . Figure M shows in fact the dielectric decrement versus its corresponding 
relaxation frequency for each dipole-active eigenmode of the given clusters. For clusters 
of two or three particles, there is one important active eigenmode close to 1/2, while for 
clusters of four particles there are two active eigenmodes. 

It can be conjectured that for a general linear cluster made of m particles, there are m — 1 
eigenvalues close to 1/2, of which the largest one is always dipole-active and has the largest 
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FIG. 5: (Color online) Effective permittivity for clusters (shown in the inset) of one, two, three, 
and four cells connected by tight and smooth junctions. The field is either parallel (solid lines 
with symbols) or perpendicular (solid lines only) to the cluster axis. The effective permittivity 
either increases strongly with the number of cells for parallel geometry, or does not change for a 
perpendicular geometry. 

weight. In fact one can show that for two cells connected by smooth and tight junction 
characterized by parameter h, (1/2 — X2) oc h 2 when h — >■ 0, or more precisely (1/2 — X2) 
is proportional with the solid angle encompassed by the missing part of a cell when it is 
connected with other cell in the dimer. The proof is based on the theorem of the solid angle 



12J. The generalization to a finite cluster is also straightforward to (1/2 — X2) oc h 2 /m (in 
that case the solid angle encompassed by the middle junction is proportional to h 2 /m. The 
weight of the second eigenmode is P%i = (x • N 1 |-u 2 )(^2| n • Ni)/^i- If we consider that 
the surface of the cluster is determined by the function g(z) then, up to a constant factor, 
(t>2|n ■ Ni) ~ g(0) 2 = h 2 for two cells connected by smooth and tight junctions. The proof 
considers that the second eigenfunction of M' is an antisymmetric combination of constant 
distributions on each part of the dimer. This assertion is confirmed in Figure [7J Moreover, 
(x-Ni|m2) jV\ is weakly dependent on h. Therefore, for a parallel setting of the field stimulus, 
P% i/(l/2 — X2)) which is the measure of the dielectric decrement of low-frequency relaxation, 
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FIG. 6: (Color online) Pk^/(l/2 — Xk) versus (1/2 — Xk) f° r clusters of up to four cells. The second 
eigenvalue has the largest contribution to intensity of relaxation. 

is finite and it increases when the number of cells is increased. The increase of relaxation 
decrement when m — > oo is physically limited by as<^S (1/2 — Xk) °\i since the membrane 
conductivity is not strictly 0. 

Due to cluster's shape and membrane properties, the variation of Pfc i/(l/2 — Xk) and 
(1/2 — Xk) with respect the eigenmode k determines a low frequency relaxation when the 
dipole-active eigenvalue Xk is close to 1/2. We note here that for dipole- active eigenmodes of 
ellip soids the term (1/2 — Xk) is called the depolarization factor and has analytical expression 
37| . Prolate spheroids with longitudinal axis much larger than the transverse axis (needles) 
have the longitudinal depolarization factor approaching and the transverse depolarization 
factor approaching 1/2. More precisely, for a long prolate spheroid, the longitudinal depo- 



larization factor scales as (1/2 — % 2 ) oc a 2 /a 2 z , a z > a a 



as a x — > 0. On the other 



hand, extensive numerical calculations support the fact that cylinders with the same aspect 
ratio behave similarly to prolate spheroids {37J1. Thus, it is not hard to observe that the 
low-frequency relaxation of linear clusters of cells connected by tight junctions is similar 
to that of a needle or a thin cylinder as long as the cluster and as thick as the junction. 
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FIG. 7: The first (dotted line) and the second (solid line) eigenfunction of M (upper panel) and 
(lower panel) for a dimer whose shape is depicted by dashed line in the lower panel. The 
second eigenfunction of is an antisymmetric combination of almost constant distributions on 
each part of the dimer. 
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On the other hand, the high-frequency relaxation of the cluster shows the relaxation of a 
suspension of spheroids with the same volume as the volume of a single cell. Therefore the 
dielectric spectrum for a suspension of clusters is the same as the spectrum of a two species 
suspension made of thin cylinders and spheroids. 



IV. CONCLUSIONS 



We present a theoretical framework based on a spectral representation of BIE and able to 
calculate the dielectric behavior of linear clusters with a wide range of shapes and dielectric 
structures. The theory agrees with the results of Pauly and Schwan for a sphere covered 
with a shell Q, IsL In fact, for spheroids, our theory is the same as the analytical results 
of Asami et al. 2]. We present extensive calculations of clusters with shapes resembling 
MDCK cells. 

A practical numerical recipe to compute the effective permittivity of linear clusters with 
arbitrary number of cells is provided. Examples are given for cluster with shapes described as 
r{6) in spherical coordinates or using (z, ip) parameters as {x = g(z) cosy?, y = q(z) sincp, z}. 



Other studies in the literature used only spherical coordinates representation 
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27, 



35|. A 



direct relation between the geometry and dielectric parameters of the cells and their dielectric 
behavior described by a Debye representation has been formulated for the first time. Other 



work [22], which is based on a closely related spectral method |18H21|. found a direct relation 



linking the geometry and electric parameters to the dielectric behavior only for homogenous 



particles. Moreover, the method used in 



22] treats only particles with spheroidal geometry. 



We show that the spectral representation provides a straightforward evaluation of the 
characteristic time constants and dielectric decrements of the relaxations induced by cell 
membrane. We prove that the effective permittivity is sensitive to the shape of the embed- 
ded particles, specially when the linear response operator has strong dipole active modes 
(with large weights pk)- A low-frequency and distinct relaxation occurs when the largest 
dipole-active eigenvalue is very close to 1/2. Clusters of living cells connected by tight 
junctions or very long cells have such an eigenvalue. Our results also shed a new light on 
the understanding of recent numerical calculations {38 ] performed with a boundary element 
method on clustered cells where the low-frequency relaxation is attributed to the tight (gap) 
junctions connecting the cells. The method used in 38| does not use the confocal geometry 
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assumption . 

The present work has several implications and applications. We emphasize the capabilities 
of dielectric spectroscopy to monitor the dynamics of cellular systems, e.g., cells during cell 



cycle division, using synchronized yeast cells 



cells 
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ay 



391 ]. or monolayers of interconnected 



4l| . Also the method is able to assess the dielectric behavior of linear aggregates 



or rouleaux of erythrocytes, where the ellipsoidal or cylindrical approximations are not 



adequate 
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43| 



The proposed representation is a powerful alternative to finite element or other purely 
numerical approaches, because it provides the analytical framework to explain and predict 
the complex dielectric spectra occurring in bioengineering applications. Extension of this 
method to other surfaces of revolution, for example linear clusters with more than 4 particles, 
is straightforward providing an adequate parametric equation is available. Finally, in many 
cases (e.g. shapes with high symmetry) the method is faster, offers accurate solutions and 
last but not least can be integrated in fitting procedures to analyze experimental spectra. 
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Appendix: Integration over <p 

The integrals over (ip — ip') are performed with the following elliptic integrals, 



T dcp 



E 



o {a-bcospf 12 y/a-ba + b \ a-b 



2b 



(A.l) 



COS ip 



o (a — b cos ip) 



3/2 



dip 



2 1 



\Ja-bb 



a + b 



-E 



2b 



K 



2b 



(A.2) 



cos 2 ip 



dip 



'o (a-bcosp) 3/2 Va-bb 2 



2a 2 - b 2 



-E 



2b 



2aK 



2b 



(A.3) 
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where K(x) and E(x) are the complete integrals of the first and second kind, respectively 
29]. 
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